Draft version February 2, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 6/22/04 



O 
O 
(N 

1—5 

o 



> 

in 

o 
o 

I 

o 



% 



ON THE PLANET AND THE DISK OF COKUTAU/4 
Alice C. Quillen, Eric G. Blackman, Adam Frank, Peggy Varniere 

Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627 
Draft version February 2, 2008 

ABSTRACT 

Spitzer observations of the young star CoKuTau/4 reveal a disk with a 10 AU hole that is most 
likely caused by a newly formed planet. Assuming that the planet opened a gap in the viscous disk, 
we estimate that the planet mass is greater than 0.1 Jupiter masses. This estimate depends on a lower 
limit to the disk viscosity derived from the time scale needed to accrete the inner disk, creating the 
now detectable hole. The planet migration time scale must at least modestly exceed the time for the 
spectrally inferred hole to clear. The proximity of the planet to the disk edge implied by our limits 
suggests that the latter is perturbed by the nearby planet and may exhibit a spiral pattern rotating 
with the planet. This pattern might be resolved with current ground based mid- infrared cameras 
and optical cameras on the Hubble Space Telescope. The required sub-Myr planet formation may 
challenge core accretion formation models. However, we find that only if the planet mass is larger 
than about 10 Jupiter masses, allowing for a high enough surface density without inducing migration, 
would formation by direct gravitational instability be possible. 

Subject headings: stars: individual (CoKuTau/4) — stars: planetary systems — planetary systems: 
protoplanetary disks — planetary systems: formation 



1. INTRODUCTION 



There remains considerable debate about the nature 
and time scales associated with planet formation in the 
disks surroundi ng young stars ( Bos^l2002: iMaver et al. I 
120021 iPollacket al. 1996) . The recent discovery of a 



age 



10'' 



years) 



with 



young star CoKuTau/4 (age t^^^ 

a disk con taining a 10 AU hole IjForrest et al.l l2?ffl^ 
iD'Alessioe t al. 2004) suggests that planet formation can 
take place quite early in the evolution of protostcUar 
systems. The critical theoretical link between the in- 
ner hole in the disk and the presence of a planet is the 
planet's ability to open a gap which impedes mass flow to 
smaller radii. Disks with edges and corresponding inner 
cleared regions have been seen in older objects on large 
scales such as HR 479 6A, with a disk edge at ^ 65 AU 
iSchneider et alJ il999l) and GM Aur with a disk edge at 
^ 800 AU. l|Rice et al.l200,'^ . However the recent Spitzer 
observations for smaller scale hole in the younger CoKu- 
Tau/4 systems provide new insights in the process which 
form planets and the evolution of planet-disk systems. 
CoKuTau/4's age and the presence of the inner disk hole 
place limits on time scales for the disk-planet interaction, 
and constrain the properties of the disk and planet. As 
we will show, the resulting estimates support the plausi- 
bility of the planet-disk scenario for CoKuTau/4. 

Because the outer disk is still present and has a spectral 
energy distribution similar to other T-Tauri stars lacking 
inner holes, we suspect that the outer disk is still accret- 
ing. Massive planet formation is unlikely to occur at ex- 
tremely small radii from the star, so following the forma- 
tion of a planet, accretion of the inner disk would leave 
a hole. Without the planet holding the accretion disk 
back, the disk edge would not be as sharp as implied b y 
the spectral energy distribution l)D'Alessio et alJl2004|) . 
For there to be an observable hole, the viscous time scale, 
Tj,, at the outer edge of the hole must be less than the 



age of the system. That the planet has not yet migrated 
to smaller radii and the presence of an inner hole allow 
us to constrain the properties of the planet and disk and 
estimate the probability of finding such systems. 

2. CLEARING THE INNER DISK AND CONSTRAINTS ON 
DISK VISCOSITY 

We parametrize the disk viscosity, v, with a Reynolds 
number, TZ = r^Q,/v where Q, is the angular rotation rate 
of a particle in a circular orbit around the star and r is 
the radius. Using a viscosity v = acgh, where a is the 
viscosity parameter, Cs is the sound speed and h is the 
vertical disk scale height, we have 
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where Vc = rf2 is the velocity of a particle in a circular 
orbit. The time scale for the disk to accrete inward is 
Ti, ^ TZTorb/'^T^ where Torb is the orbital period. To es- 
timate Tjy we require the disk aspect ratio h/r. From 
hydrostatic equilibrium, we have h/r ~ Cs/vc. The 
sound speed is given by c^ ~ kBT/rUp where ks is Boltz- 
mann's constant, and rup is the proton mass. The tem- 
perature T, if set by radiative balance with the star, is 

/ ^ \i/4 
T ^ le 2 , where a is the Stefan-Boltzmann con- 

\ iDTrr^cr y ' 

slant. Heat released in the disk by accretion could lead 
to a higher temperature than that estimated above. Al- 
ternatively, the disk could be self-shielding, cooling the 
mid-plane below the temp erature estimated above (e.g., 
iChiang fc Goldreichlll997|) . Keeping these compHcations 
in mind we use this temperature as a starting point. 

CoKuTau/4 is an Ml. 5 star at a distance of 140pc 
with an estimated mass of O. SMfD, Tgne. = 10^ 
years, and luminosity of .6^0 f ljD'Alessio et alJ 120041: 
iKenvon fc HartmannI Il995^ ) . Assuming an inner edge 



of 10 AU and the luminosity and mass given above. 

Electronic address: [aquillen, afrank, blackman, pvarni ©pas.rochester.edaF, the disk edge we find T -- 80K, Cs ~ 0.8km/s and 



Vc ^ 7kni/s. The orbital period at 10 AU is ^ 40 years. 
Using the this sound speed and hydrostatic equihbrium, 
h/r ~ 0.1. However, if the disk is optically thick, h/r 
could be lower. We use h/r = 0.05 below, though this 
particular choice is not significant. 

If a planet exists in CoKuTau/4's disk at r =10 AU, 
then the hole formed because disk material is prevented 
from accreting across the planet's orbital radius by the 
transfer of orbital angular momentum from the planet to 
the disk. The presence of the hole therefore implies that 
the material in the inner disk (r <10 AU) has had time 
to accrete onto the star. This requires t^ < Tage — IMyr, 
which in turn implies 

r. = TlT^rbl'^^ = a"' (J^) rort/2^ < IMyr. (2) 

Using Torb =40 years, the above inequality implies 
that TZ < 1.6 X 10^. Using h/r — 0.1, we also re- 
quire that a > 6 X 10~^. These values are within the- 
oretical expectations. A viscosity parameter a ~ 0.01 
typically emerges from protostellar disk model fitting 
(jHartmann et al.lll998t) . Note that if a or /i were orders 
of magnitude lower, the inner disk would not have had 
time to accrete onto the star. In principle, the viscosity 
of the inner disk could differ from that of the outer disk, 
or either could have dropped in the past million years, 
but here we are just assuming the simplest steady case. 
The result highlights that the hypothesis that the inner 
disk has accreted onto the central star is consistent with 
age constraints and reasonable disk properties. If instead 
we use a k, 0.01 and h/r = 0.05 we find 7?. ~ 4 x lO** and 
Tjy w 2 X lO^yr, also less than Tage and thus in principle 
consistent with the observations. However, the planet 
formation time must be less than the viscous time at the 
radius where the planet forms. A viscous time signifi- 
cantly less than 1 Myr challenges core accretion models 
if the planet formed at the corresponding radius. 

Instead of accreting onto the central star, the inner disk 
could have been depleted by multiple planet formation, 
or agglomeration into large dust grains and planetesimals 
within 10 AU. Here however we focus on the simplest 
paradigm of a single planet and viscous accretion. 

3. GAP OPENING AND CONSTRAINTS ON PLANET MASS 

The inner disk will begin accreting onto the star once 
a newly formed planet opens up a gap. The condi- 
tion for opening a gap provides a limit on the planet's 
mass. Without a gap, the disk would accrete unimpeded 
through the orbital radius of the planet. A gap decouples 
the inner and outer disk, except through their interaction 
with the planet. Because the planet acts as a time de- 
pendent gravitational potential perturbation, it can res- 
onantly drive waves at Lindblad resonances into a the 
gaseous outer planetesimal d isk ( Goldrcich & Trcmaine 
Il978t iLin &: PapaloizoullTOTQ : Ward 1997). These waves 
carry angular momentum and th erefore govern both 
how a planet opens gaps (e.g., iBrvden et al.l 119991: 
lArtvmowicz fc Lubow"1994') as well as the the radial mi- 
gration (e.g., Nelson et al. 2000; Ward 1997). 

To open a gap, a planet must be sufhciently massive 
that spiral density waves dissipa ted in the disk overcome 
the inward fiqw due to viscosity ( ijLin fc Papaloizoull97g: 
IBrvden et a l."1999':'Ward [ll997|) '). This leads to the con- 
dition l|Nelson et a l. 20(^ 

q>407e"\ (3) 



where q = AIp/AI^, the mass of the planet divided by 
that of the star. For 7?, < 1.6 x 10^ estimated above 
for CoKuTau/4, the gap opening condition implies that 
g>2.5xl0~'' 
Jupiter. 



or Mp > 0.1 Mj where Mj is the mass of 



4. MIGRATION, MASS LOSS AND SURFACE DENSITY 

Once formed, a planet's interaction with a surround- 
ing disk may lead to the transfer of angular momentum 
and t he migration of the planet closer to the star l|Wardl 
Il997|) . For more massive planets, migration can occur 
after a gap is opened (denoted Type II migration). The 
condition for opening a gap and the nature of Type II mi- 
gration are linked. A gap is maintained when the torque 
density from spiral waves driven at different resonances 
balances the inward torque due to viscous accretion. If 
the planet mass is less than or comparable to the disk 
mass with which it interacts then the planet will migrate 
on viscous time scale; it behaves a s "just anoth er par- 



ticle" in disk and r, 



mig 



(e.g., INelson et all 12000). 
However this result is subject to the surface density pro- 
file; outwa rd migration rather than i nward migration 
may occur (|Masset fc Papaloizoull200^ . 

If the planet mass is large compared with the disk mass 
with which it interacts then the inertia of the planet slows 
migration and Tmig > t^- That CoKuTau/4's planet still 
resides a large distance from the star implies that sig- 
nificant migration has yet to occur. The lack of signifi- 
cant migration along with the presence of a hole implies 
Tciear < Tmig, whcrc Tdear is the timcscalc for the gas 
within the planet's semi-major axis to accrete onto the 
star. If we use Tdear ^ t^, then we require Tmig > ^V for 
type II migration. 

After a gap is opened, the inner disk accretes onto the 
star while outer disk experiences a pile up of material at 
the edge of the disk exterior to the planet. Once M^dge ~ 
Mp, migration would take place on a viscous time scale. 
Here M^dge = T^Eg where Eg is the surface density just 
outside the disk edge. The accretion rate of the outer 
disk may be crudely estimated to be 
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(4) 
Significantly larger accretion rates would have led to 
larger accumulations of mass in the disk edge and an ear- 
lier onset of migration. We note that accretion rates in 
the range of value of Ma = 10~^ — 10"^" Mq/ji are con- 
sistent with observations of accretion rates in million year 
old evolved T-Tauri systems (jCalvet et al.ll200ry . Larger 
values of Ma could be accommodated in our calculations 
by taking a larger planet mass planet and/or a larger re- 
quired mass in the disk edge to initiate migration. If we 
assume that the build up of M^dge > Mp would have led 
to inward migration, its absence constrains the disk sur- 
face density through Alp > irr'l'Ee, where rp is the semi- 
major axis of the planet. Using our limit Mp > 0.1 Mj 
estimated above, we find Eg < 4 gm cm^'^ in the disk 
edge, a plausible value for disks around young stars. 

The discovery of a planet at ^ lOAU from CoKuTau/4 
is therefore consistent with the hypothesis that inward 
disk migration has not proceeded to completion because 
sufficient mass has not accumulated in the disk edge. 



5. PREDICTED DISK MORPHOLOGY 

Having established that a young sub-Jovian mass 
planet orbiting at 10 AU in the CoKuTau/4 disk is plau- 
sible, we now consider the gravitational and hydrody- 
namic interaction between the plan et and di s k. Us i ng the 
hydrodynamics code developed by ' Masseti l)2000l 12002(1 
we have performed a 2D hydro simulation using the pa- 
rameters estimated in the previous sections. Figure 1 
shows the morphology from a simulation with planet 
mass ratio q — 3 x 10^^, Reynolds number TZ = 10''', 
and disk aspect ratio h/r = 0.05. In the simulation, 
the planet was initially set into a circular orbit, with the 
disk edge located at 1.1 times the semi- major axis of the 
planet's orbit. The initial surface density was taken to 
be E = So(l.lrp/r)-i where Eq = IQ-^M^/r^ is the 
surface density at the disk edge and Vp is the planet's 
semi-major axis. For r < l.lvp we set the disk density to 
be 100 times lower than that at the disk edge to approx- 
imate an initial inner hole. The planet is free to migrate 
via gravitationally induced angular momentum exchange 
with the disk, and can accrete gas within its Roche Lobe. 
Note that the simulated disk is not massive enough for 
self-gravity to play a role. Figure 1 shows the gas density 
at time t — lOOrort. after the beginning of the simulation. 

The principle conclusion from our simulations relates 
to spiral density waves driven into the disk from inter- 
actions with the planet. Because of the proximity of the 
planet to the disk edge, the disk contains more than one 
Lindblad resonance. Multiple spiral density waves can be 
driven in the disk edge by the planet at these resonances. 
The range in radius wh ere the waves are la unched is of or- 
der the scale height, h l|Artvmowic2ll99,'^ . Consequently 
we expect the winding of the spiral pattern to depend on 
the scale height and hence on the disk temperature, i.e. 
the spiral wave would be more tightly wound if the disk 
were cooler and the scale height smaller. The two-armed 
pattern rotates with the planet and is probably caused by 
the combined effect of more than one density wave, as ex- 
plained by ( Qgilvie fc Lubow 2002) . The spiral pattern 
could be observable by high angular resolution imaging. 
At a distance of 140pc, 0.1" ~ 14 AU. The features ex- 
hibited in the simulation are then close to the resolution 
limit of ground based 10-meter class telescope at 10/im 
or by the Hubble Space Telescope (HST) at optical wave- 
lengths. Thus it is possible that the spiral pattern might 
also be detected in optical scattered light images. In 
addition, the asymmetry of the disk hole caused by the 
arm which extends toward the planet may be detectable 
through imaging. 

When the Reynolds number or the planet mass is 
higher, the disk edge would be further away from the 
planet. Then the lower m Lindblad Resonances are the 
dominant sites of density wave driving. Consequently, 
2 or 3 density peaks, corresponding to the constructive 
interference of 2 and 3 armed waves might be detected 
in images of the disk edge. If the disk aspect ratio is 
higher than considered here, a one-armed pattern rotat- 
ing with the planet dominates over the two-armed one 
seen in Figure 1. 




Fig. 1- Morphology of a planet disk system with 
planet mass of 0.1 Jupiter mass for the CoKuTau/4 
system. The disk aspect ratio is h/r = 0.05 and 
Reynolds number TZ = 10^. Because of the proxim- 
ity of the planet and disk, high order spiral density 
waves are driven by planet planet at Lindblad reso- 
nances. If these waves constructively interfere, the 
result is a spiral pattern that rotates with the planet. 
The planet in this simulation is accreting gas. 

Hydrodynamic simulatio ns have shown differences be- 
tween 2D and 3D models l|Makita et al.l l2'600'). particu- 
larly in the opening angle of the spiral density waves. In 
a 3D disk, the waves may not constructively interfere as 
they do in our 2D simulati on, becau se of differences in 
the dispersion relation (Qgilvie fc Lubow 2002) . Further 
3-D simulations are required to better determine struc- 
tures in the disk formed via planet-disk interactions. 

In the simulation shown, we allow the planet to accrete 
80% of the material found within its Roche Lobe. The 
remaining 20% can flow p ast the planet into the inner 
disk (as seen previously bv lLubow et al.lll9 9&). account- 
ing for the non-zero gas density within the planet's semi- 
major axis. If the disk has the high viscosity considered 
here, then it is likely that the planet is still accreting 
significantly. A planet that is still accreting may be sur- 
rounded by a hot observable circumplanetary disk (e.g., 
iLubowet anil999() . 

The Spitzer observation s imply that there is v ery little 
dust within the disk edge IjD'Alessio et al."2004"l. though 
the presence of larger bodies is less constrained. Better 
hydrodynamic modeling will help to understand the flow 
past the planet into the inner hole and consequently on 
the disk, planet mass and planetary accretion. 

6. DISCUSSION AND CONCLUSIONS 



The recent discovery by Forrest et al. ( 2004) of a young 
system with a 10 AU hole (..D'Alessio et al...2004) allows 
us to explore new constraints on the evolution of young 
planets and circumstellar disks. Based on the assumption 
that the inner disk has accreted onto the star CoKuTau/4 
within a time equivalent to the age of the star we esti- 
mated that TZ < 1.5 X 10^. Using this, we evaluated the 
condition for a planet to open a gap in the disk and found 



that the planet within the disk edge of CoKuTau/4 is 
greater than O.lMj. The apparent lack of inward migra- 
tion of the planet leads to estimates of the disk accretion 
rate and surface density which are consistent with obser- 
vations of evolved T-Tauri systems. Given the inferred 
planet mass, we expect the disk edge to be very near the 
planet. This implies that the planet could be accreting 
material and would interact strongly with the disk by 
driving waves into the disk from resonances. Idealized 
simulations suggest that this may produce a two-armed 
spiral pattern, rotating with the planet and extending a 
few disk scale heights away from the radius of the planet. 
Structure in the disk could be detectable with 0.1" high 
resolution imaging in scattered optical light with HST or 
by ground based 10m class telescopes in the mid-infrared. 
The contemporaneous presence of both a sharp edge 
and an inner hole (suggestive of viscous inner disk clear- 
ing) implies that the hole clearing timescale is less than 
the planet migration timescale, Tdear < ^mig, otherwise, 
the planet would have disappeared along with the hole, 
and the hole re-filled. The probability to see both an 
edge from the planet induced gap AND an disk hole for 
r < {lOAU - Srgap/2) is then 



P = [T„,^g{10AU) - T,{10AU - Srgap/2)]/Tag 



(5) 



where Svgap is the gap width opened initially by the 
planet and the time scales are meant to be taken at 
the radius in parentheses. For Svgap « 10 AC/, P ^ 
(TrmgilOAU) - Tyis{10AU))/Tage- Ouce the planet ap- 
pears at lOAU, it takes at least one viscous timescale 
for material to pile at the disk edge, and presumably 1 
more for migration to take place, so then P ^ T,j/Tage- 
For a = 0.01 the estimate for t^ of section 2 would then 
give P ~ 10%. This argument presumes that planet ap- 
pearance at IQAU occurs before a viscous time at IQAU . 
(This is more restrictive than the minimal condition that 
a planet must form faster than a viscous time scale at its 
formation radius because the planet could have formed 
at r > lOAC/.) 

More robust processes to explain why the inward mi- 
gration ti me scale Tm.in < Tricar are pos sible and may be 
necessary. IMasset fc Papaloizoul (J2003J) find that surface 
density profiles shallower than r~^/^ can induce an out- 
ward planet migration rather than inward which could 
certainly explain the presence of wall and disk. Alterna- 



tively, if the inner accretion disk had a radially depen- 
dent viscosity coefficient, or incurred a change of accre- 
tion mode (e.g. non-viscous transport) the surface den- 
sity profile might evolve so as to produce a hole within 
the planet radius even if the planet migrated on a local 
viscous time. More work on these possibilities for CoKu- 
Tau/4 will be needed. 

Finally, we note that the planet-disk scenario for 
CoKuTau/4 implies that planets of greater than two 
Neptune masses can form in < 10^ years. Th is may 
push the envelope of core-accretion models (iPoUack et alJ 
I1996J) and this may appear to favor pl anet f orma- 
tion by direct gra v itatio nal instability l|Bossl I2002L 
120031 iMaver et all l2002|) . Our results constrain 
the properties of the disk from which the planet 
formed: The Toomre instability parameter for a ther- 
mally supported disk at the planet orbit radius Vp is 

Q - 200 (^) {j^y\j^y'{!^), where 
we have scaled to to the values based on our lower mass 
limit for the planet in CoKuTau/4. Planet formation 
by direct gravitational instability is only possible for 
Q < 1.5. For this condition to be satisfied, the disk 
would have had be remarkably thin and/or have a very 
high E in the not-to distant past. Our upper limit on 
S is proportional to Mp, so that only if Mp = lOAfj, 
would the upper limit of E < 400g/cm^, allow Q to be 
low enough for planet formation by direct gravitational 
instability. Observational constraints on the surface den- 
sity are needed. 
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